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Abstract 

We have measured the branching fractions of the colour suppressed decays B'^ — > 'ip(2S)K*~^ (892) 
and tfj (2 S) K*^ (892) using a data sample of 84 million BB events recorded by the Belle 

detector on the T(46') resonance. The branching fractions for the charged and neutral mode 
decays are (8.13 ± 0.77 ± 0.89) x 10"^ and (7.20 ± 0.43 ± 0.65) x lO"'', respectively. 

PACS numbers: 13.25.Hw, 14.40.Nd 



INTRODUCTION 

The decays B V^(2S') K* are colour-suppressed two-body decays with two vector mesons 
in the final state. Measurements of such processes can provide insight into the interplay 
between the weak and strong interactions. 

We study eight different decay channels, which are combinations of two decay modes 
for the ^(2^): ij{2S) U and ij{2S) J/^jti+ti-; and four different K* final states: 
K*~^ — ^ K+vr", K*^ Ks7!"^, K*^ K^n^, and K*^ Ksn^ (charged conjugate processes 
are implied throughout this paper). The branching fractions for B ip{2S) K* have been 
measured previously by the CLE0j3| and CDPj^ experiments. Our much larger data sample 
permits us to make more precise determinations of these quantities. 

DATA SAMPLE AND EVENT SELECTION 

This measurement uses a data sample of approximately 84 million BB events, corre- 
sponding to an integrated luminosity of 78 fb~^, collected at the T(4S) resonance by the 
Belle detector jsj at KEKB[J]. KEKB is an energy- asymmetric double storage ring that 
collides 8 GeV electrons with 3.5 GeV positrons. 

The Belle detector is a general-purpose large-solid- angle spectrometer with a 1.5-T 
solenoidal magnet. It is designed to study the properties of B mesons, charmed hadrons, 
and T leptons, and also to study two-photon interactions. Charged particles are detected 
by a silicon vertex detector (SVD) that surrounds the interaction point, and then tracked 
by a cylindrical wire drift chamber (CDC). The CDC provides a three-dimensional momen- 
tum measurement of each track. Charged particles are identified via information from silica 
aerogel Cherenkov counters (ACC), a time-of-fiight scintillator barrel (TOF), and specific 
ionization {dE/dx) measurements in the CDC. Electromagnetic showers are detected in an 
array of CsI(T£) crystals (ECL) that is located just inside the solenoidal coil. Muons and 
Kl mesons are identified by glass-electrode resistive-plate counters (KLM) embedded in the 
solenoid's magnetic fiux return. 

Hadronic events are selected from the raw data sample by requiring at least three good 
charged tracks emerging from near the e~^e~ interaction point, within a cylindrical region 
of r < 2.0 cm and \z\ < 5.0 cm, aligned along the positron beam axis, and at least one 
cluster of 0.1 GeV or more within the ECL barrel. In the center-of-mass (CM) frame, the 
component of the total momentum along the beam axis of all charged tracks and neutral 
ECL clusters must be below 0.5y/s/c, where y/s = 10.58 GeV is the total available energy. 
In the laboratory frame, the visible energy (the sum of the charged track momenta and 
the neutral cluster energies) must exceed 0.2^/s and the neutral-only part must be between 
O.ly^ and 0.8y/s. Continuum (non-resonant e+e" — >• qq) events are suppressed relative to 
e'^e^ T[4:S) by requiring that the ratio R2 of the second to the zeroth Fox- Wolfram 
momentsjBI falls below 0.5. 

Candidate charmonium mesons are reconstructed through the decays J/ip i~^£~, 
ip{2S) — >■ and ipi^S) — > J/ipTT^n' i~^i~n~^7r~, where i represents either an electron 

or a muon. Electrons are identified by a maximum-likelihood technique that uses energy 
and momentum information from the ECL and CDC, respectively. If there is a neutral 
ECL cluster within 0.05 rad of the electron/positron direction, it is assumed to arise from 
bremsstrahlung, and its energy is added to that of the lepton. Both tracks are required 
to originate from a common vertex by requiring that their distance of closest approach, 



measured along the beam axis, satisfies \dz\ < 5.0 cm. Candidate J/ip and 'ip{2S) mesons 
are retained if tlie mass difference between tlie reconstructed e^e~ invariant mass and the 
nominal J/^p or ip{2S) mass falls in the range -0.150 < {M^+e- - Mj/^) < 0.036 GeV/c^ 
or —0.186 < (Mg+e- — M(/i(2S)) < 0.036 GeV/c^. (The windows are wider on the low side 
of the charmonium mass peak to allow for residual bremsstrahlung.) Muons are identified 
by a maximum likelihood technique that uses the momentum information from the CDC 
in combination with the range and transverse scattering information from the KLM. As 
for electrons, muon pairs must satisfy the common-vertex requirement of \dz\ < 5.0 cm. 
Candidate J/ip and ip{2S) mesons are retained if the reconstructed fi~^fj^~ mass satisfies 
-0.060 < (M^+^- -Mj/^) < 0.036 GeV/c^ or -0.070 < (M^+^- - M^(2S)) < 0.036 GeV/c^. 
The invariant- mass spectrum of the candidate inclusive ip{2S) sample is plotted in Fig. ^ 




FIG. 1: Invariant mass spectra of ip{2S) candidates. The left plot shows ip{2S) — > while right 
plot shows ip{2S) J/ipir^TT^. 

Additional 4'{2S) candidates are reconstructed through their decay to J/ipTi^ii~. The 
J/ip meson is reconstructed via its decay to dileptons, as described above. Here, the mo- 
menta of the leptons and the parent J/ ip candidate are redetermined in a mass-constrained 
fit assuming the nominal J/ip meson mass. Pions are selected from charged tracks after 
electron and muon rejection. To reduce combinatorial background, the vr+vr" invariant mass 
must exceed 0.4 GeV/c^, and the CM momentum of the ip{2S) candidate must fall below 
1.7 GeV/c. The three-dimensional distance between the reconstructed vertices of the lepton 
pair and the pion pair must be less than 0.5 cm. We retain a candidate ip{2S) meson if its 
mass exceeds the reconstructed J/ip mass by (0.59 ± 0.01) GeV/c^. The spectrum of the 
mass difference between the ip{2S) candidate and the J/ip daughter is shown in Fig. [T] 

We reconstruct K*^ candidates from their decays to K~^7i~ or Kgir'^, and K*^ candidates 
from their decays to K~^tt^ or KsTt"^. Charged tracks are identified as kaons if they satisfy 
the kaon selection criterion derived from a maximum likelihood technique that distinguishes 
them from pions using information from the CDC, ACC, and the TOF; otherwise, they 



are classified as pions. Neutral pions are reconstructed from pairs of photons (i.e., neutral 
ECL clusters) whose individual energies exceeds 60 MeV and whose invariant mass is in 
the window 0.12 < M^^ < 0.15 GeV/c^. (The photon and pion momenta of the surviving 
candidates are also recalculated assuming the nominal mass of the 7r° meson.) 

Ks mesons are reconstructed from two oppositely charged pions whose distance of clos- 
est approach to the interaction point exceeds 0.25 mm in the plane perpendicular to the 
beam axis. If the two tracks have SVD hits, they must originate from a common ver- 
tex that is within the SVD but detached from the interaction point. Otherwise, the 
azimuthal coordinate of the common vertex must agree with the azimuthal direction of 
the pion pair within 0.1 rad. The invariant mass of the pion pair must be in the range 
0.482 GeV/c2 < M^^ < 0.515 GeV/c^. 

When combining a neutral pion with either a charged or neutral kaon to form a K* 
candidate, we require the K* helicity angle (the angle of the kaon with respect to the K* 
direction in the K* rest frame) to satisfy cos 6* < 0.85. This requirement provides rejection 
of the slow-pion background. The invariant mass of the K^7t~ or Kstt~^ combinations 
and the Kstc^ or K^tt^ combinations must agree with the nominal K* meson mass within 
0.075 GeV/c^ and 0.090 GeV/c^, respectively. For the K*^ K^it~ mode, the momenta of 
the kaon and pion are recalculated in a vertex-constrained fit. 

MEASUREMENT OF BRANCHING FRACTIONS 

Candidate B mesons are reconstructed by combining il^{2S) candidates with K* candi- 
dates. For each B candidate, we calculate the momentum ps and energy Eb in the CM 
frame, and from these the beam-constrained mass Mbc = \l -E'^eam ~ Pb ^^cl the energy dif- 
ference Ai? = Eb — -Ebeam- (-Ebeam = -\/^/2 is the beam energy in the CM frame.) 

If there are multiple B candidates in an event, the best one is chosen by first discarding 
all but the best-reconstructed iIj{2S) candidate, and then selecting the B candidate with 
the smallest ■ This quantity is defined by = Y^iij^i — f^iY/'^iy where fii is the central 
value of the measured parameter rrii — i.e., the mass of the reconstructed ip{2S) or K* — 
and CTj the corresponding uncertainty. The B candidate selected in this way is retained 
only if it lies within a specified region of the Mbc-AE plane: —0.2 < AE < 0.2 GeV and 
5.2 < Mbc < 5.3 GeV/c2. 

A signal box, \AE\ < 0.03 GeV (0.04 GeV for tt^ modes) and 5.27 < Mbc < 5.29 GeV/c^, 
was determined using a large sample of Monte Carlo signal events. We have verified that the 
Mbc and AE resolutions agree well between data and simulation. In the cases where the final 
state contains only charged tracks [i'^i'K^TT" and i~^i~7T'^TT~ K~^tt~), the track momenta 
are recalculated assuming a common vertex constraint; this improves the resolutions in Mbc 
and AE for signal B^ candidates while having no effect on the background distribution. 

Using the Mbc and AE signal box cut to restrict the sample of B meson candidates, 
we obtain Kn mass distributions for K*'^ and K*^ candidates. These distributions are 
fitted to Breit-Wigner p-wave functions for signal and second-order polynomial functions for 
background. The Kn mass spectra of the exclusive K* candidates are shown in Fig. |21 -B — >■ 
ip{2S)K* backgrounds are classified as coming from the qq continuum, combinatorics, feed 
across, and non-resonant Kit decay. The relative importance of each background category 
varies with the signal mode under consideration. Backgrounds of each type were studied for 
each mode using a combination of Monte Carlo samples and data sidebands, as described 




FIG. 2: Invariant mass of the (a) K*^ candidates reconstructed from K vr^ and Ksir^ pairs, (b) 
K*^ candidates reconstructed from K^tt^ and Kst^"^ pairs. 



below. 

We simulate the continuum background using the proportion 1 : 1.6 : 1 for b : uds : c 
quarks. A Monte Carlo sample with a size corresponding to twice the data sample was 
generated and subjected to the selection criteria outlined above. Since no events survive our 
cuts, we conclude that the continuum contribution is negligible. 

Combinatorial backgrounds are analyzed by using the sidebands of the experimental data. 
(A MC study was also considered, but owing to complications arising from non-resonant B — »• 
iI){2S)Kt: backgrounds and contributions from higher K* resonances, it was not adopted.) 
We take the areas in the AE-Mbc plane given by 5.2 < M^c < 5.27 GeV/c^ and -0.2 < 
AE < -0.05 GeV and 0.05 < AE < 0.12 GeV for our sideband analysis. We exclude 
the 0.12 < AE < 0.2 GeV region from this study to avoid possible non-combinatorial 
contributions from ip{2S)K and ip{2S)n decays. Background from feed across is studied 
by using exclusive signal MC data generated using the program EvtGen j^. The generated 
particles are tracked through the detector using the program GEANT Eight 20,000-event 
samples, corresponding to the eight B — > ip{2S)K* modes studied, are used. For each decay 
mode studied, feed across backgrounds from the remaining seven modes are summed. The 
total amount of feed across relative to signal yields is 8.1-12.9%, depending on the mode. 

Finally, backgrounds from non-resonant B iI){2S)Ktt modes are studied by using a 
sideband in the Ktt mass spectrum with masses in the range 1.1 < Mk-k < 1-3 GeV/c^. We 
obtain an estimated background contribution from non-resonant Kn modes of about 6.8% 
for V(25)K*+ and 4.7% for ij{2S)K*'^. 

Measurements of the branching fractions are based on fits to the beam-constrained mass 
(Mfjc) distributions. Separate fits are performed for each of the eight submodes under study. 
These distributions comprise signal and background events as discussed in the previous sec- 
tions. In the fits, the signal is described by a Gaussian whose mean and width are determined 
from MC. The combinatorial background is described using an ARGUS function whose 
shape parameters are fixed using sideband data. The feed across background is modeled 
by a Crystal Ball function where both the shape and the level are fixed by MC. The 
level of the feed across background depends on the B — >■ ip{2S)K* branching fractions being 
measured, so an iterative procedure is adopted wherein the branching ratios obtained using 
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FIG. 3: Beam-constrained mass plots and fits for — > tp{2S)K*^ . The left-hand plots are 
for ip{2S) — > and the right-hand plots are for ip{2S) J / ipir^ ir^ . The top plots show 

K*~^ — > K~^Tr^ and the bottom plots show K*~^ — > KsTr~^- The dashed lines show the combinatorial 
contribution to the background. The dotted lines (green) show the feed across background and 
the solid lines with hatched areas (red) show the contributions from non-resonant decay. The gray 
lines give total background. 



initial estimates for the feed across modes are used to carry out a second analysis. Since the 
feed across backgrounds are a small fraction of the signal, this procedure quickly converges. 
The non-resonant Kir background is described by a Gaussian of fixed shape and amplitude. 
The level of the combinatorial background is allowed to vary. Figs. El and |3] show the fit 
results for — > iIj{2S)K*~^ and —>■ %1){2S)K*^ decays, respectively. Fitting results from 
the Mfoc distributions are checked by carrying out a second set of fits to the AE distributions. 
This method yielded consistent results, with yield differences varying from 1 to 3%. The re- 
construction efficiency for each decay mode is obtained using signal MC data with the same 
cut criteria. Table H] shows the number of signal events, the reconstruction efficiency and 
the total number of background events for each decay mode. World average values |lO| are 
used for the branching fractions for ip{2S) ipi^S) — > J/ip ti'^ti' , J/ip — >• and 

K* —* Kit decays. Systematic errors arise from uncertainties in detector efficiency, particle 
identification efficiency, reconstruction efficiency, and beam-constrained mass fitting. We es- 
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FIG. 4: Beam-constrained mass plots and fits for — > ■tp{2S)K*^. The left-hand plots are 
for '4>{2S) — > and the right-hand plots are for ip{2S) J/ipTT~^TT~. The top plots show 

K*^ K'^TT~ and the bottom plots show K*^ — > Kstt^. The dashed lines show the combinatorial 
contribution to the background. The dotted lines (green) show the feed across background and 
the solid lines with hatched areas (red) show the contributions from non-resonant decay. The gray 
lines give total background. 



timate these as follows: tracking efficiency (1.0% per track); particle identification efficiency 
(3.0% for leptons and less than 1% for kaons and pions); background contribution (5.0 - 7.0% 
decay mode dependent); neutral pion reconstruction (4.0%); number of BB events (1.0%); 
branching fraction for secondary decay (1.7%); and the effect of polarization (1.7%). The 
uncertainty in the background, which is the dominant contribution to the overall systematic 
error, is due to uncertainties surrounding the contributions of higher-mass K* resonances to 
the 1.1 < Mktt < 1.3 GeV/c^ mass region used to estimate the non-resonant contributions 
to the Kit mass region under the i^'*(890) peak. 

The measured branching fractions are summarized in Table HI Table ITTl shows the branch- 
ing ratios that result from the weighted average of the four final states for the charged 
and neutral mode decays. Table |n] also shows previous measurements by CLEO Q for 
comparison. 



Decay channel 


M2S) 


sig 


Sr-fo ( /O ) 


iV hhn 


B X 10"^ 


ip{2S)K*+{K+7r'^) 


a 


50.6 ±8.6 


15.85 


33.2 


7.96 ±1.41 ±0.84 




tpTTTT 


47.6 ± 9.8 


5.25 


57.1 


8.97 ± 1.78 ±0.96 


i){2S)K*+{KsTT+) 


U 


32.2 ±6.0 


15.5 


9.8 


7.69 ±1.46 ±0.82 




■IpTTTT 


36.1 ± 7.0 


6.42 


18.2 


8.21 ±1.66 ±1.12 


iP{2S)K*^{K+TT-) 


U 


181.4 ± 15.2 


29.5 


62.6 


7.66 ±0.62 ±0.61 




tpTTTT 


164.2 ± 14.6 


12.1 


91.6 


6.70 ±0.65 ±0.74 


il){2S)K*^{KsTT^) 


U 


13.2 ±4.2 


10.3 


7.3 


9.35 ±2.73 ±0.97 




■IpTTTT 


7.9 ±4.3 


4.01 


13.4 


5.75 ± 3.43 ± 1.02 



TABLE I: Number of signal events, reconstruction efficiency, number of background events and 
measured branching fractions for each decay mode. 





Belle (B x 10^) 


CLEO {B X 10^) 


Decay channel 


B+ ^{2S)K*+ 


8.13 ±0.77 ±0.89 


8.9 ±2.4 ±1.2 




7.20 ± 0.43 ± 0.65 


7.6 ± 1.1 ± 1.0 



TABLE II: The results of measured branching fractions after taking weighted average for different 
decay channels. 

CONCLUSION 

The branching fraction for B — > iIj{2S)K* (890) has been measured using a sample of 84 
million B pairs. The results are consistent with earlier measurements from CLEO f]\ and 
CDF 1^], but offer improved statistical accuracy. 
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